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We investigate rogue wave solutions in a three-component coupled nonlinear Schrodinger equation. 

With the certain requirements on the backgrounds of components, we construct a new multi-rogue 
wave solution that exhibits a structure like a four-petaled flower in temporal-spatial distribution, 
in contrast to the eye-shaped structure in one-component or two-component systems. The results 
could be of interest in such diverse fields as Bose-Einstein condensate, nonlinear fibers and super 
fluid. 
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I. INTRODUCTION 

Rogue wave(RW) is localized in both space and time, 
which seems to appear from nowhere and disappear with¬ 
out a trace M- It is one of the most fascinating phe¬ 
nomena in nature and has been observed recently in non¬ 
linear optics [3] and water wave tank 4]. The studies of 
RW in single-component system have indicated that the 
rational solution of the nonlinear Schrodinger equation 
^NLS) can be used to describe the phenomenon well [5|— 

A variety of complex systems, such as Bose-Einstein 
condensates, nonlinear optical fibers, etc., usually involve 
more than one component [9]. Recent studies are ex¬ 
tended to RWs in two-component systems (Ml. Some 
new structures such as dark RW have been presented 
numerically 0 and analytically [l2j |. Moreover, it was 
found that two RWs can emerge in the coupled system, 
which are quite distinct from high-order RW in one- 
component system [l2|. In the two-component coupled 
systems, the interaction between RW and other nonlin¬ 
ear waves is also a hot topic of great interest dlil El. 
It was shown that RW attracts dark-bright wave in |9[. 

In the present paper, we further extend to investigate a 
three-component coupled system considering the number 
of the modes coupled in the complex systems is usually 
more than two. With the certain requirements on the 
backgrounds of components, we construct a new rational 
solution that can be used to describe the dynamics of 
one RW, two RWs, and three RWs in the system. One 
structure like a four-petaled flower is found in the coupled 
system: there are two humps and two valleys around 
a center in the temporal-spatial distribution, which is 
quite distinct from well-known eye-shaped one presented 
before. We discuss the possibilities to observe them in 
nonlinear fibers. 


The paper is organized as follows. In Section II, we 
present exact vector RW solutions and the explicit con¬ 
ditions under which they could exist. The dynamics of 
them are discussed in detail. In Section III, the possibil¬ 
ities to observe them are discussed. The conclusion and 
discussion are made in Section IV. 


II. EXACT VECTOR ROGUE WAVE 
SOLUTIONS AND THEIR DYNAMICS 

It is well known that coupled NLS equations are often 
used to describe the interaction among the modes in non¬ 
linear optics, components in BEC, etc.. We begin with 
the well-known three-component coupled NLS, which can 
be written as 

i ~W + l)x^ + 41V’ 1 ! 2 + hfol 2 +1^31 2 ] = o, 

j ‘~tt + ~dx ^ + 41V’ 1 ! 2 + + \1/J3\ 2 ]lp2 = o, 

*77 + tSt + 4IV’i| 2 + \H 2 + = o, (i) 

where ipj (j=1,2,3) represent the wave envelopes, t is 

the evolution variable, and x is a second independent 
variable. The Eq.(l) has been solved to get vector soli- 
ton solution on trivial background through Horita bilin¬ 
ear method in 0. Performing Darboux-transformation 
from a trivial seed soltuion with ^3 = 0, one could get the 
bright-bright solitons in 0- It has been reported that 
solitons could collide inelastically and there are shape¬ 
changing collisions for coupled system, which are differ¬ 
ent from uncoupled system 0. However, it is not possi¬ 
ble to study vector RW on trivial background. Next, we 
will solve it from nontrivial seed solutions. The nontrivial 
seed solutions are derived as follows 

Vho = si exp [i 2 (s\ + s\ + s\)t + ik\x — ik\t\, (2) 

^20 = s 2 exp [i2{s\ + si + s\)t + ik 2 x — ik\t\, (3) 

^30 = S 3 exp [i2(s\ + si + s\)t + ik%x — ik^t] (4) 
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where Sj(j=l,2,3) are arbitrary real constants, and de¬ 
note the backgrounds in which localized nonlinear waves 
emerge, Aq, k 2 and k% denote the wave vectors of the 
plane wave background in the three components respec¬ 
tively. From physical viewpoint, the relative wave vector 
value is important. One of the three components can be 
seen as a reference to define the wave vectors of the other 
two. Then, we can set &2 = 0 without losing generality. 
To derive the rational solutions, we find that there are 
some requirements on the amplitude of each component, 
and the difference of their wave vectors should be related 
with the amplitude in a certain way. The conditions un¬ 
der which one can get vector RW with no other type 
waves are 



, , , V2 

ki = -k 3 ,k 3 = —si, 

V2 

S 2 = -^-Sl,S3 = Si- 


(5) 

(6) 


With the conditions and s 1 = 1, the generic form of 
vector RWs could be given as 



Vh = 

(•- 

^2 = 

(*- 

^3 = 

(*- 

The expressions 

for j 

part. It is seen 1 

that 

the expressions, 

A A 

ters. Therefore, 

the 1 


Hi(x,t) 

Gi(x,t) 

H 2 (x,t) 

G 2 (x,t) 

H 3 (x,t) 

G 3 (x,t) 


r .9 1 . x n 

exp \i - 

PL 2 V2 J 
exp [5 it] 

T/I - ’ 

r .9 1 . X 

exp b-+*- 7 =] 


(7) 

(8) 
(9) 


RWs solution, which can be verified by the following RWs 
plots. There are mainly three cases for the generalized 
vector RW solution. 

One vector rouge wave — When As = 0, A* = 0, only 
one vector RW can be observed in all components, as 
shown in Fig. 1. Interestingly, we find that there is a 
novel shape for vector RW solution. The density distribu¬ 
tion shapes of the localized waves in ^1 and ^3 are quite 
different from the well-known eye-shaped one. There are 
two humps and two valleys around a center, and the cen¬ 
ter’s value is almost equal to that of the background, 
as shown in Fig. 1(a) and (c). This structure can be 
called as four-petaled structure in temporal-spatial dis¬ 
tribution. Moreover, the humps or valleys in 'ipi corre¬ 
spond to the valleys or humps in ^ 3 . However, the den¬ 
sity distribution in ^2 is similar to the eye shaped RW in 
single-component system, as shown in Fig. 1(b). There¬ 
fore, the whole density is still the well known eyes shape, 
as shown in Fig. 1(d). The novel shape should come from 
the cross phase modulation effects, since the shape can 
not be observed in scalar RW. For two-component cou¬ 
pled systems, it has been found numerically that there 
are dark RWs in one component of the coupled system in 
El- The dark RW has been given exactly in our previ¬ 
ous paper in fl2j |. Based on these results, we expect that 


FIG. 1: (color online) The evolution plot of one RW in cou¬ 
pled system, (a) for one RW in component, (b) for one RW 
in component, (c) for one RW in ips component, (d) for the 
whole density of the three components. It is seen that the den¬ 
sity distribution shapes of RW in ipi and ^3 are different from 
the eyes shape. The coefficients are A\ = 0, A = 1, A 3 = 0, 
and A 4 = 0. The ”L” and ”S” in color bar denote large and 
small values of density respectively. This holds for all figures. 


there should be some novel structures in more than three 
modes coupled systems. 

Two vector rouge waves — When A 4 = 0, there are 
two vector RWs appearing in the temporal-spatial distri¬ 
bution, shown in Fig. 2 . When A\ < 0, there are two 
vector RWs emerging at a certain propagation distance, 
as shown in Fig. 2(a)-(c). It is seen that the structures of 
the two RWs in every mode are similar, and only the sizes 
are different. The four-petaled structure RWs emerge in 
Vq ,3 and the eye-shaped ones emerge in ^2 component. 
When Ai » 0, the two distinct RWs emerge at different 
propagation distances, shown in Fig. 2(d)-(f). There is 
a rotation on x — t distribution plane for the two RWs. 
Varying the parameter A, 3 , one can observe the interac¬ 
tions between the two RWs. 

Three vector rouge waves — When A 4 7 ^ 0, there are 
three vector RWs appearing in the temporal-spatial dis¬ 
tribution, as shown in Fig. 3. When A 4 A 2 < 0 and 
| H 2 1 >> | Ai|, there are three vector RWs emerging very 
clearly at a certain propagation distance, as shown in 
Fig. 3(a)-(c). The structures of them in each compo¬ 
nent are similar with different sizes. The valleys in A 
still correspond to the humps in i/j 3 component. There 
are three eye- shaped RWs in A- When A 4 A 2 > 0 and 
IA 2 I >> | Ai|, the three RW emerge clearly at different 
propagation distances, as shown in Fig. 3(d)-(f). The 
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FIG. 2 : (color online) The evolution plot of two vector RWs in coupled system, (a) for the RWs in ipi component, (b) for the 
RWs in ip 2 component, (c) for the RWs in ^3 component. The coefficients are A 1 = 0, A 2 = 40, A 3 — 8 , and A 4 = 0. (d), (e) 
and (f) show the evolution of RWs in 'ip 1 , ^ 2 , and ip 3 respectively. The coefficients are A± = 150, A 2 = 40, A 3 = 8 , and A 4 = 0. 


RW seems to have changing velocity on retarded time. 
Varying the parameters, we can investigate the interac¬ 
tion between these vector RWs conveniently. For exam¬ 
ple, making A 2 approach A 4 , we can observe the inter¬ 
action of the three RWs, such as Fig. 4. It shows that 
RW’s shape can be changed greatly through interacting 
with others. Two RWs almost fuse to one valley in ^3 
component with the condition, as shown in Fig. 4(c). 

It should be noted that the distribution shapes of the 
two and three RWs in the whole temporal-spatial dis¬ 
tribution are very distinct from the high-order RW in 
one-component system presented in [3, [RJ, [lfj. In the 
one-component systems, it is not possible to observe just 
two RWs appearing in the whole temporal-spatial distri¬ 
bution even for high-order RWs. Three RWs can emerge 
on temporal-spatial distribution for second-order scalar 
RW, but their distribution shapes are different from the 
three RWs obtained here. Therefore, we call the gener¬ 
alized vector RW solution as multi-vector RW. 


III. POSSIBILITIES TO OBSERVE THESE 
VECTOR NONLINEAR WAVES 

Considering the experiments on RW in nonlinear fibers 
@|, which have shown that the simple scalar NLS could 


describe nonlinear waves in nonlinear fibers well, we ex¬ 
pect that these different vector RWs could be observed in 
three-mode nonlinear fibers. One could introduce three¬ 
mode optical signals into a nonlinear fiber operating in 
the anomalous group velocity dispersion regime, marked 
by Vb U = 1,2,3). Firstly, we assume the nonlinear 
coefficients for the three modes are equal. The disper¬ 
sion and nonlinear coefficients are 1 and 2 in normalized 
units respectively. Then, the backgrounds can be given 
by Eq.(5) and ( 6 ). Explicitly, the amplitude of ipi is set 
to be s 1 , the amplitude of ^2 is <$2 = -^si, and the am¬ 
plitude of ^3 is 53 = si. The wave vector of ^2 marked 
by &2 is set to be zero, and the wave vectors of ipi and 
'ips are k\ = —£ 3 , and ks = ^s\. Then the initial opti¬ 
cal signals are given by the presented vector RW solution 
Eq.(7)-(9). For nonlinear fibers, the coordinates x and 
t should be changed to be t and z, which are retarded 
time and propagation distance respectively. Under the 
corresponding conditions, one, two, or three vector RWs 
could be observed in the nonlinear fiber. 

On the other hand, it is well known that there are 
spatial optical solitons in planar waveguide. The sim¬ 
ilar conditions can be derived directly through coordi¬ 
nates transformation for the coupled NLS in a planar 
waveguide E3- The RWs in multi-mode planar waveg- 




































4 



-20 -10 


10 20 




-20 


-10 


0 

X 


10 20 


(a) 


(b) 


(c) 


15 


10 


-10 


-is [ 

-10 -50 5 10 

X 




0 5 10 

X 


(d) 


(e) 


(0 


FIG. 3: (color online) The evolution plot of three vector RWs in coupled system, (a) for the RWs in t/u component, (b) for the 
RWs in ip 2 component, (c) for the RWs in ^3 component. The coefficients are A 1 = 0, A 2 = 20, A 3 = 0, and A 4 = — 1. (d), (e) 
and (f) show the evolution of RWs in t/u, ^ 2 , and ^3 respectively. The coefficients are A± = 0, A2 = 20, A3 = 0, and A4 = 1 . 


uide could be observed too. Additionally, considering the 
studies on multi-components Bose-Einstein condensates 
EU 23 , we expect that the vector RWs could also be 
observed with the condition. 


IV. DISCUSSION AND CONCLUSION 

In summary, we find some novel structures for vector 
RWs in a three-component coupled system. The novel 
shape lies in that there are two humps and two valleys 
around one center in the temporal-spatial distribution. 
The obtained multi-vector RW solution can be used to 
describe one, two, and even three vector RWs in the cou¬ 
pled system, whose density distribution shapes are differ¬ 
ent from high-order scalar RW. The corresponding condi¬ 
tions for their emergence are presented explicitly. Under 
these conditions, the ideal initial signals for them can be 
given from the generalized rational solution, which could 
be helpful for experimental observation. Based on these 


results, we expect that abundant novel structures could 
exist in more-than-three-component coupled systems. 

The coupled system can be used to describe three- 
component BEC, multi-mode optical transmission, and 
so on. Therefore, we believe that these nonlinear waves 
can be observed experimentally. As an example, the pos¬ 
sible way to observe vector RWs in three-mode nonlinear 
fiber is discussed here. Recently, the high-order modu¬ 
lation instability with RWs has been observed in non¬ 
linear fiber optics in [21]. There are many possibilities 
to observe similar phenomena in coupled nonlinear fiber 
system. 
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Appendix A: The analytic expressions for Hj(x,t) 
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Hi = —(1 + *) ^ 6 iAi + (6 — 6*)j4 3 — 6 i\/ 2 A 4 + (6 + 12 i)A 3 t + 6 \Z 2 A A t — 6 iA 3 t 2 — (6 — 6 i)V 2 A A t 2 + 2 s/ 2 A 4 t 3 
+A 2 ( 6\/2 — 6v^2 1 + 6 ix) + 6 V 2 A 3 X + (6 — 6*)j4 4 x — 6 \p 2 A 3 tx + (6 + 12 i)A 4 tx — 6 iAH 2 x + 3 iA 3 x 2 
+ 3 V 2 A 4 X 2 — 3\/2A 4 tx 2 + *^4 4 x 3 J |Vj V2A 1 + (6 + 12i)V2A 3 t — (6 + 6*)j4 4 i — 6\/2A 3 t 2 + (18 — 6*)j4 4 i 2 
+4iAit 3 — (6 — 6 *)j 43X — 12*4l3te + (6 + 12 i)V 2 A A tx — 6 V 2 A A t 2 x + 3V2A 3 x 2 — (3 — 3*)j4 4 x 2 — 6 iA A tx 2 
+V 2 A 4 X 3 + 6^4 2 (—1 + i — 2 it + a/ 2 x)J , 

Gi = V 2 ^36Al + 362I3 + 54v / 2^ 3 ^4 + 54 + 2 + 108Alt 2 - 36\f2A 3 A A t 2 - 72^3^ + 48V2A 3 A 4 t 3 - 24A 2 A t 3 

+36Alt 4 - 48v / 2^4 3 A i t£ 4 + 108^^ 4 - 24A^ 5 + 8 Alt 6 + 36V2A%x + 108A 3 A 4 x + 54\/2A 2 x - 72\[2A%t 2 x 
+144A 3 A A t 2 x - 36V2Alt 2 x - 48A 3 A 4 t 3 x + 48V2A\t 3 x + 24A 3 A 4 t 4 x - 48V2A\t 4 x + 36Alx 2 + 54 ^ 2+3 Aix 2 
+54^4 4 x 2 — 36j4 3 fx 2 + 36A 3 t 2 x 2 — 72V2A 3 A 4 t 2 x 2 + 72^4 4 i 2 x 2 — 24^4 2 i 3 x 2 + 12^4 4 i 4 x 2 + 36j4 3j 44X 3 
+18V2A 2 4 x 3 - 24A 3 A4tx 3 + 24A 3 A A t 2 x 3 - 24V2A 2 A t 2 x 3 + 9A 2 x 4 + 9A 2 x 4 - 6 A^x 4 + 6 A 2 4 t 2 x 4 + 6 A 3 A 4 X 5 
— 12+ 2 (— 3A4 + 6A 4 t 2 — 8A 3 t 3 + 4 ^ 4 ^ — 3A 4 \/2x — 3A 4 x 2 — A4X 4 + 6A 3 V2t 2 — 6A 3 t 2 x + 6A 3 tx — 3V%A 3 ) 
+ 12 j 4 i( 6 j 42 X + 6A 3 t — 6A 3 t 2 + 3 ^ 3 X 2 + 6\Z2A 4 t 2 — 6A 4 t 2 x + 6 A 4 ix + 2 I 4 X 3 ) 

+72A 2 A 3 x + 36A24I3X 3 + A^x 6 + 36A 2 (1 - 2 1 + 2 t 2 + x 2 )] . 


H 2 = 72V2A 2 + 72iV2A% + (72 - 72i)V2A!A 3 + 144 iA 2 A 3 + (144 - 144*)Ai^ 4 + I44iy/2A 2 A4 - (144 + 144*) i/2A%t 

+(144 + \44i)V2A 1 A 3 t + (216 + 72i)V2A 2 3 t - (288 + 144*) V2A 2 A 4 t + (288 + 144 i)A 3 A 4 t - 144 V2A%t + 144V2Alt 2 
-144V2A 1 A 3 t 2 - 288A 2 A 3 t 2 + 216i^2A 2 3 t 2 + 288A 1 A 4 t 2 + (144 - 288i)V2A 2 A 4 t 2 - 144^ 3 ^4< 2 
+(144 - 360i)V2A 2 A t 2 - (144 + 144*) V2A%t 3 + (192 + 192i)V2A 2 A 4 t 3 + (192 + 192 i)A 3 A 4 t 3 + (192 + 48*)^2^4t 3 
+72V2A 2 3 t 4 - 96V2A 2 A 4 t 4 - 192^ 3 ^4i 4 + (96 + 120i)V2Alt 4 - (48 + 48i)V2A 2 4 t 5 + 1 6V2A 2 A t 6 + 144V/2A1212X 
+(72 + 72*)v / 2 j 4 2 j4 3 x + 144*+ 2 x + (72 - 72i)V2A 1 A a x + 144A 2 A a x + 144 *a/ 2213214X - (144 + 144*) 

+(144 + 144*)ViL4i7Ux + 144V2A 3 A A tx + (288 + 144 i)A 2 A tx + 144V2^ 2 ^ 3 t 2 x - 288 A%t 2 x 
— 144v / 2^4ij44< 2 x + (144 + 144*)\/2j43^44£ 2 x — 144^4£ 2 x — (96 + 96*)v / 2^4 3 j44t 3 x + (192 + 192*)^44£ 3 x 
+48 v / 2j4 3 j44< 4 x — 192^4 2 £ 4 x + 72\Z2^4|x 2 + 72 a/2j4ij4 3 x 2 + (36 + 36i)V2A 3 x 2 + 72\/2j42^44X 2 
+(72 + 144*)+3^4X 2 + 72iV2A 2 4 x 2 - (72 + 72i)V2A 2 3 tx 2 + 72^ A 2 A tx 2 + 72s/2A 2 3 t 2 x 2 + 12V2A 2 A t 2 x 4 
—288A 3 A 4 t 2 x 2 + (72 + 72i)V2A 2 A t 2 x 2 - (48 + 48i)V2A 2 A t 3 x 2 + 24V2A 2 A t 4 x 2 + 72\[2A 2 A 3 x 3 + 2 V 2 A 2 A x 6 
+24V/2A1414X 3 + (48 + 24*) V2+3-44X 3 + (24 + 48*)24 2 x 3 - (48 + 48*)^44 3 Aite 3 + 48V2A 3 A A t 2 x 3 
—96 A 4 t 2 x 3 + 18\Z2 j4 2 x 4 + 24\/2^42^44X 4 + (12 + 6*)%/2^4 2 x 4 — (12 + 12i)s/2A 4 tx 4 + 12v / 2j4 3j 44X 5 , 

G 2 = 36V2Af + 36^2^1 + 72^ 2 ^s + 36V2A 2 3 + 36V2^ 2 ^ 4 + 108^3^4 + 54^2^1 - 72^2^|t 
+ 72 V 2 .A 1 A 3 t + 72V2Alt 2 - 72V2AiA 3 t 2 - 144 A 2 A 3 t 2 + 108V2A 2 3 t 2 + 12\Z2A 2 A a x 4 
+144^i^4< 2 - 72\[2A 2 A A t 2 - 72A 3 A 4 t 2 - 72\f2A%t 3 + 96\f2A 2 A A t 3 + 96 A 3 A A t 3 - 24^2 A 2 A t 3 
+36\/2 A 2 3 t 4 - 48\Z2A 2 A 4 t 4 - 96A s A 4 t 4 + 108^2 A\t 4 - 24^2 A 2 A t 5 + 8 V2A 2 4 t 6 + 72T2^i^ 2 x 
+72\/2 j42j 43X + 72^4 3 x + 72A 2 A a x + 108v / 2^43j 44X + 108^4 2 x — 72 v / 2 j 42j43<x + 72\/2^4i^44ix 
+72\/2 j42j43< 2 x — 144^4 3 < 2 x — 72\/2^4 ij44< 2 x + 144\/2j43j44£ 2 x _ 72^4 2 < 2 x — 48"\/2j43j44i 3 x 
+96 A\t 3 x + 24\f2A 3 A A t 4 x - 96A 2 A t 4 x + 36\Z2+ 2 x 2 + 36 a/)L4i413X 2 + 36V2A3X 2 + \l2A\x 6 
+36v / 2^4 2 Aix 2 + IO84I3+4X 2 + 54V2A 2 a x 2 - 36V2A 2 3 tx 2 + 36V2A 2 3 t 2 x 2 - 1442l32l4t 2 x 2 
+72v / 2A 2 t 2 x 2 - 24^L4 2 fV + 12V24l 2 t 4 x 2 + 36s/2A 2 A 3 x 3 + 12 a/ 2 +iAix 3 
+36v / 2A 3 A 4 x 3 + 36 A 4 X 3 - 24v / 2 J 43+4te 3 + 24s/2A 3 A A t 2 x 3 - 48A 2 4 t 2 x 3 + 9V2+ 2 x 4 
+9V2A 2 x 4 - 6\/2A\tx 4 + 6V2Alt 2 x 4 + 6V2A 3 A a x 5 . 
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H 3 = (1 + i) j 6 A — 6 iA 3 t — 6 A 3 t 2 + 6 V 2 A 4 I 2 — 2 iV 2 A 4 t 3 + 6 iV 2 A 3 tx — QiA^tx — 6 Aii 2 x + 3 A 3 X 2 + 3 iV 2 A 4 tx 2 
+A 4 X 3 + 6 A 2 (iV 2 t + x^jj ^6%/2 j4i + (6 + 12 i)V 2 A 3 t — (6 + 6 i)A 4 t — 6 \Z 2 A 3 t 2 + (18 — 6i)Ai 2 + 4*Aii 3 
— (6 — 6 *)j43X — 12 *j 431 x + (6 + 12 i)V 2 A 4 tx — §\p2A/tf? x + 3V2A 3 x 2 — (3 — 3*)Ax 2 
— QiA 4 tx "I - \I 2 A 4 X 64 I 2 (—1 4~ 1 ) — 2 it \/ 2 x^ j , 

G 3 = V2 [36,41 + 36^| + 54 V 2 A 3 A 4 + 54 A 2 + 108,4|l 2 - 36V2A 3 A 4 t 2 - 72A 2 3 t 3 

+48V2A 3 A 4 t 3 - 24A\t 3 + 36A 2 3 t 4 - 48V2A 3 A 4 t 4 + 108^ 2 ! 4 - 24,4 2 1 5 + 8 A\t & + 36V2A 2 3 x 

+108 A Ax + 54\/2 A\x - 72V2A 2 3 t 2 x + 144AA* 2 x - 36\/2 A\t 2 x - 48AAi 3 x 

+48y/2A4t 3 x + 24+ 3 Al 4 x — 48\Z2^4 2 £ 4 x + 36,4 3 x 2 + 54y / 2AAx 2 + 54,4 2 x 2 — 36A 3 tx 2 

+36,4 3 i 2 x 2 - 72V2A 3 A 4 t 2 x 2 + 72 A 2 4 t 2 x 2 - 24A 2 4 t 3 x 2 + 12,4li 4 x 2 + 36A 3 A 4 x 3 

-12A (A(-3 + 61 2 - 8 1 3 + 41 4 - 3V2.X - 3x 2 - x 4 ) + A 3 (6t 2 V2 - 6t 2 x + 6 tx - 3V2 - 6x - 3x 3 )) 

+12A ^6Ax + A 3 (6t — 6t 2 + 3x 2 ) + A(6i 2 \/2 — 6t 2 x + 6 tx + x 3 )^ 

+ 18V2Alx 3 - 24A 3 A A tx 3 + 24A A*V - 24^2,4 2 i 2 x 3 + 9^ 3 x 4 + 9,4 2 x 4 - 6A 2 lx 4 
+6^1 2 x 4 + 6 A 3 A 4 X 5 + Alx e + 36Al(l -21+ 2 1 2 + x 2 )] . 



